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SUMMARY 

Relations  -were  analytically  obtained  for  the  prediction  of  radial 
distributions  of  velocity  and  temperature  for  fully  developed  laminar 
flow  of  gases  and  of  liquid  metals  in  tubes  with  fluid  properties 
variable  along  the  radius.  The  relations  are  applicable  to  both  heat¬ 
ing  and  cooling  of  the  fluid.  By  use  of  the  relations  for  velocity  and 
temperature  distributions ,  relations  were  obtained  among  Nusselt  number , 
friction  parameter,  and  ratio  of  wall  to  bulk  temperature. 

The  Nusselt  number  and  friction  parameter  were  found  to  be  inde¬ 
pendent  of  Reynolds  number  and  Prandtl  number.  The  effects  of  ratio  of 
wall  to  bulk  temperature  on  Nusselt  number  and  friction  parameter  could 
be  eliminated  by  evaluating  the  fluid  properties  at  specified  tempera¬ 
tures  in  the  fluid. 


INTRODUCTION 

Most  of  the  analyses  of  laminar  flow  and  heat  transfer  in  tubes, 
for  instance  the  uniform -wall -temperature  solution  given  in  reference  1, 
have  been  carried  out  under  the  assumption  that  the  fluid  properties  do 
not  vary  with  temperature.  Somewhat  more  work  on  variable  fluid  prop¬ 
erties  has  been  done  for  the  laminar  boundary  layer  of  a  flat  plate. 

(See,  for  instance,  references  2  and  3).  In  reference  4,  velocity 
distributions  for  laminar  flow  in  tubes  with  variable  viscosity  and  some 
approximate  relations  for  temperature  distributions  and  heat-transfer 
coefficients  are  given. 

In  the  analytical  investigation  reported  herein,  which  was  con¬ 
ducted  in  the  NACA  Lewis  laboratory,  both  velocity  and  temperature  dis¬ 
tributions  are  obtained  for  fully  developed  laminar  flow  in  tubes  of 
gases  and  of  liquid  metals  with  variable  fluid  properties  .  These 
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relations  are  applicable  to  both  heating  and  cooling  of  the  fluid.  By 
use  of  these  velocity  and  temperature  distributions,  Nusselt  numbers  - 
and  friction  parameters  are  obtained.  The  results  should  be  applicable 
to  flow  at  points  where  the  ratio  of  distance  from  entrance  to  diameter 
of  the  tube  is  on  the  order  of  one— twentieth  of  the  Reynolds  number  or 
greater. 


D 

d 

f 

*0 

8 

h 

k 

kb 

k0 

m 

Nu-b 

m0 

Pr 


SYMBOLS 

The  following  symbols  are  used  in  the  report: 

specific  heat  of  fluid  at  constant  pressure,  (Btu)/(rb)(°F) 
inside  diameter  of  tube,  (ft) 

exponent,  value  of  which  depends  on  variation  of  viscosity  of 
fluid  with  temperature 

friction  factor,  2^^/ (pu-^) 

friction  factor  with  density  evaluated  at  tx 

friction  factor  with  density  evaluated  at  wall  temperature 

acceleration  due  to  gravity,  (32.2  ft/sec2) 

heat-transfer  coefficient,  -  t-^),  (Btu) /(sec) (sq  ft)(  F) 

thermal  conductivity  of  fluid,  (Btu) (ft) /(sec) (sq  ft)(°F) 

thermal  conductivity  of  fluid  evaluated  at  t^, 

(Btu) (ft)/ (sec) (sq  ft)(°F) 

thermal  conductivity  of  fluid  evaluated  at  tQ, 

(Btu) (ft)/(sec)(sq  ft)(°F) 

Nusselt  number,  hD/k 

Nusselt  number  with  conductivity  evaluated  at  t-^ 

Nusselt  number  with  conductivity  evaluated  at  tx 
Nusselt  number  with  conductivity  evaluated  at  tQ 
Prandtl  number,  Cpgp /k 
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q.  rate  of  heat  transfer  toward  tube  center  per  unit  area, 

(Btu)/(sec)(sq.  ft)) 

rate  of  heat  transfer  at  wall  toward  tube  center  per  unit  area, 
(Btu) / (sec) (sq  ft)) 

gf  dimensionless  heat -transfer  parameter,  q_/q.Q 

r  radius,  distance  from  tube  center,  (ft) 

Tq  inside  tube  radius,  (ft) 

rT  dimensionless  radius  parameter,  t/tq 

Be  Beynolds  number,  pu^D/p 

Bex  Beynolds  number  with  density  and  viscosity  evaluated  at  tx 

BeQ  Beynolds  number  with  density  and  viscosity  evaluated  at  wall 

t  absolute  temperature,  (°B) 

t-^  bulk  or  average  static  temperature  of  fluid  at  cross  section 

of  tube,  (°B) 

tx  temperature  in  fluid,  x(tQ  -  t-^)  +  t^,  (°E) 

tQ  absolute  wall  temperature,  (aR) 

tf  dimensionless  temperature  parameter,  t/tg 

t-^1  dimensionless  bulk-temperature  parameter,  t-^/tg 

u  velocity  parallel  to  axis  of  tube,  (ft/sec) 

u-^  bulk  or  average  velocity  at  cross  section  of  tube,  (ft/sec) 

u r  dimensionless  velocity  parameter,  pQu/(TQr0) 

u^1  dimensionless  bulk  velocity  parameter,  ^oub^rOrO^ 

x  number 

z  distance  from  tube  entrance,  (ft) 

\  dimensionless  parameter,  kQt0/(q.oro) 
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dimensionless  friction  parameter ,  T oro/ (ut)P ) 

dimensionless  friction  parameter  with  viscosity  evaluated  at  tQ 
dimensionless  friction  parameter  with  viscosity  evaluated  at  tx 
absolute  viscosity  of  fluid,  (lb-sec/sq  ft) 

absolute  viscosity  of  fluid  evaluated  at  t-^,  (lb-sec/sq.  ft) 
absolute  viscosity  of  fluid  at  wall,  (lb-sec/sq  ft) 
mass  density,  (lb -sec2 /ft4) 

bulk  or  average  density  at  cross  section  of  tube,  (lb -sec2 /ft4) 

mass  density  of  fluid  at  wall,  (lb-sec2/sq  ft4) 

shear  stress  in  fluid,  (lb/sq  ft) 

shear  stress  in  fluid  at  wall,  (lb/sq  ft) 


ANALYSIS 

In  the  present  analysis,  equations  for  the  velocity  and  temperature 
distributions  in  a  laminar  stream  of  gases  or  of  liquid  metals  through 
a  tube  are  derived  for  the  fully  developed  boundary  layer.  From  these 
equations,  relations  for  the  Nusselt  number  and  friction  factor  are 
obtained . 


Flaw  of  Gases 


The  analysis  for  gases  applies  when  the  following 


Assumptions . 
assumptions  hold: 

1  The  velocity  u  at  any  given  distance  from  the  wall  is  inde¬ 
pendent  of  distance  along  the  tube.  This  assumption  holds  for  fully 
developed  flow  when  the  rate  of  density  change  in  the  direction  along 
the  tube  is  snail.  The  effect  of  density  change  becomes  less >“*£*£* 
as  the  Telocity  Is  increased.  Laminar  flew  can  exist  at  these  yeloci 

if  the  tube  diameter  is  small. 

2.  The  temperature  difference  tQ  -  t  at  any  distance  from  the 
wall  is  independent  of  distance  along  the  tube.  This  con  l  ion 

”  .  1  •  .1 _ •  1 _ U,'  -TV 


fS  fSlyl”d  temperature  distentions  with  uniform  longitudinal 
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distribution  of  heat  transfer  in  the  tube  when  the  fluid  properties  are 
constant  and  holds  for  variable  properties  when  the  heat -transfer  rate 
varies  slightly  along  the  tube. 

3.  Frictional  heating  effects  can  be  neglected.  These  effects  are 
small  when  the  Mach  number  for  the  flow  is  small  compared  with  unity; 
the  Mach  number  is  generally  small  when  laminar  flow  through  a  tube 
exists . 


4  4.  The  viscosity,  the  thermal  conductivity,  and  the  density  vary 

with  temperature;  the  specific  heat  and  the  Prandtl  number  are  constant. 
The  variations  with  temperature  of  the  specific  heat  and  the  Prandtl 
number  of  gases  are  of  a  lower  order  of  magnitude  than  the  variation  of 
viscosity,  thermal  conductivity,  and  density. 

5.  Heat  transfer  by  conduction  in  the  gas  in  the  direction  along 
the  tube  can  be  neglected.  This  heat  transfer  is  generally  small  com¬ 
pared  with  that  transferred  by  convection  along  the  tube  except  for  very 
small  velocities. 

Velocity  distributions .  -  For  laminar  flow  the  shear  stress  is  pro¬ 
portional  to  the  velocity  gradient,  or 


The  relation  for  the  variation  of  shear  stress  with  radius  for  unaccel¬ 
erated  flow  (assumption  l)  is  obtained  by  equating  the  shear  forces  to 
the  pressure  forces  acting  on  a  cylinder  of  fluid  of  arbitrary  radius 
and  differential  length  (reference  5).  This  relation  gives 


T 


(2) 


From  equations  (l)  and  (2) 


r 

r0 


T0  ~  =  “M-  ~ 


du 

dr 


(3) 


In  order  to  convert  this  equation  to  dimensionless  form,  the  follow- 
ing  quantities  are  introduced: 


u*  s 


'cFo 


(4) 


r'  = 


r0 


(5) 
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Then 


du'  M’O  du 
drT  ~  t  0  dr 


Equation  (3)  becomes ,  when  written  in  dimensionless  form, 


The  variation  of  the  viscosity  of  a  gas  or  liquid  metal  with 
temperature  can  he  shown  from  viscosity  data  to  he  given  hy  an  equation 
of  the  form 


^0  VO 


Substituting  this  value  for  ju/m-q  in  equation  (7)  and  separating 
variables  results  in 

r  * 

du'  = - y  dr' 

t'd 

Integration  of  this  equation  between  a  point  in  the  fluid  and  the  wall 
results  in 


In  order  to  solve  this  equation,  the  relation  between  f  and  r'  must 
be  obtained.  This  relation  is  obtained  in  the  following  section. 

Temperature  distributions  .  -  For  obtaining  temperature ^distri-^ 
buttons,  a  heat  balance  is  first  written  for  an  annulus  of  fluid  wi 
inside  radius  r,  outside  radius  r  +  dr,  and  length  dz.  The  heat 
entering  the  outer  cylinder  of  the  annulus  is  2it(r  +  dr)(q  +  dq)dz  and 
Se  Sit  Saving  the  inner  cylinder  is  2„  r  g  in.  net  heat  trans- 

ferred  into  the  annulus  through  the  cylinders  is,  if  differentials  of 
higher  order  than  the  first  are  neglected,  2it  dz(rdq  +  qdr)  or 
2rt fS(rq) /Srldr  dz.  The  heat  transferred  through  the  ends  of  the  annulus 
by  conduction  is  neglected  (assumption  5) .  The  heat  picked  up  by  fluid 
passing  through  the  annulus  is  2jt  r  dr  pgu  cp  or 
2«  r  dr  dz  pg  u  cp  dt/dz.  The  heat  balance  can  then  be  written  as 
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S(rq)  St 

sT =  rpg  u  °p  ST  (10) 

The  heat  "balance  for  the  whole  cross  section  of  the  tube  is 
2n  rQ  q0  =  it  rQ2  Ph  g  ub  cp  Stb/dz 

or 

2%  =  r0  S  uh  cp  Stb/Sz  (11) 

From  assumption  2, 

S(tp  -  t)  =  o 

Sz 


or 


5tp  at  5tb 
Bz  dz  5z 


(12) 


Dividing  equation  (10)  by  equation  (ll)  and  using  condition  (12)  and 
constant  Cp  (assumption  4)  give 


1  S(r'  q/qp) 

2  ^rt 


_P_  iL 
Pb  ub 


(13) 


The  static  pressure  is  uniform  across  the  tube.  From  the  perfect  gas 
law, 

(14) 

Pb  t 

Substituting  equation  (14)  in  equation  (13 )  and  integrating  produce 
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Another  expression  for  q'  is  obtained  as  follows;  The  law  for 
heat  conduction  is 

,  dt  r- 

q  =  k  —  \- 

u  dr 

Equation  (16)  can  he  written  in  dimensionless  form  as 

,  k  dt' 

q-'  “  k  k0  dr'  ^ 


where 


^0^0 


With  the  assumptions  of  constant  Erandtl  number  and  constant  specific 
heat  (assumption  4),  the  law  for  the  variation  with  temperature  of  the 
conductivity  must  be  the  same  as  the  law  for  variation  with  temperature 
of  the  viscosity,  or 

=  ( i-Y  =  t,d  (19 


k0  VO 


Substitution  of  equation  (19)  into  equation  (17)  gives 

^  dtT 

q'  = 

Equating  equation  (15)  to  equation  (20)  gives 


_1_  _  _2_ 
r T  r 1 


liip'. xfi  £1 

t  T  U’  1  n  ¥ 1 


This  equation  can  "be  integrated  to  give 


1  /u  -  /  Eli' dr '1  dr 


1'V  °\/r.  >  vJr 


In  order  to  solve  equation  (23.),  tt,',  up,',  and  u'  must  be 

obtained  in  terms  of  t'  and  r' . 
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r1 

u-u'  =  2  /  u'r'  dr' 

u  Jo 

From  equations  (23)  and  (25) 


( 25) 


(26) 


When  equation  (26)  and  the  value  for  u' 
tuted  in  equation  (21) ,  there  results 


from  equation  (9)  are  substi- 
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Equation  (27)  gives  the  relation  between  t'  and  r'  for  various 
values  of  \  for  laminar  flow  of  gases.  The  equation  can  he  solved  by 
the  method  of  iteration;  that  is,  an  assumed  relation  between  t'  and 
r1  is  substituted  into  the  right-hand  side  of  the  equation  and  new 
values  are  calculated  by  solving  the  equation.  The  new  values  are  then 
substituted  into  the  equation  and  the  process  is  repeated.  After  the 
relation  between  t'  and  r'  for  each  value  of  \  has  been  obtained 
from  equation  (27),  the  relation  between  u'  and  r'  can  be  obtained 
from  equation  (9) . 


Flow  of  Liquid  Metals 

Assumptions.  -  Assumptions  1,  2,  3,  and  5,  which  were  used  in  the 
analysis  of  laminar  flow  of  gases,  are  also  used  in  the  analysis  of 
laminar  flow  of  liquid  metals  in  tubes.  Assumption  4  is  altered,  how¬ 
ever,  inasmuch  as  the  variations  with  temperature  of  the  density  and  the 
thermal  conductivity  for  most  liquid  metals  are  small  compared  with  the 
variations  of  viscosity.  In  the  present  analysis,  it  is  assumed  that 
only  the  viscosity  varies  with  temperature. 

Yelocity  distributions.  -  In  the  derivation  of  equation  (9)  for 
velocity  distributions  in  gas  flow,  variable  density  and  thermal  con¬ 
ductivity  did  not  enter  the  analysis.  Equation  (9)  therefore  applies 
also  to  liquid  metals . 

Temperature  distributions.  -  The  analysis  of  temperature  distribu¬ 
tions  for  liquid  metals  is  the  same  as  tnat  for  gases  except  that  P/Pb 
in  equation  (13)  and  k/k0  in  equation  (17)  are  both  set  equal  to  1. 

The  equation  for  temperature  distributions  is,  with  these  alterations, 
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This  equation  can  be  solved  in  the  same  wa y  as  equation  (27).  The 
dimensionless  hulk  temperature  t^T  for  liquid  metals  is  calculated 
from 


rather  than  from  equation  (23 ),  which  was  used  for  gases. 


(29) 


Nusselt  Numbers  and  Friction  Parameters 

Nusselt  numbers.  -  The  Nusselt  number  with  the  conductivity  evalu¬ 
ated  at  the  wall  temperature  is  defined  as 


where 


-  qo 

h  =  t0  - 


(30) 


(31) 


From  equations  (30)  and  (31)  and  the  definition  of  X.  in  equation  (18), 
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The  conventional  friction  factor  with  density  evaluated  at  the  wall 
fQ  can  "be  shown,  from  the  definitions  of  fg,  5g,  and  Reg,  to  he 

given  hy 


f 


(40) 


or 


(41) 


Conventional  Velocity-  and  Temperature -Distribution  Parameters 

The  conventional  velocity -distribution  parameter  u/u-^  can  be 
written  as 


u  _  u1 

ub  V 


(42) 


where  u'  is  calculated  from  equation  (9)  and  u^'  is  calculated  from 
equations  (25),  (9),  and  (27)  for  gases,  or  equations  (25),  (9),  and  (28) 
for  liquid  metals . 


Similarly  the  temperature -distribution,  parameter  (tg-t) /(tg-t^)  can 
be  written  as 


t0~  t 
t0"tb 


1-V 


where  t'  is  calculated  from  equation  (27)  for  gases  or  equation  (28) 
for  liquid  metals,  and  t^*  is  calculated  from  equation  (23)  or  (29). 


Special  Case  of  Uniform  Fluid  Properties 

Equations  (9)  and  (27)  or  (28)  can  be  integrated  for  the  special 
case  where  the  fluid  properties  are  constant  by  setting  the  value  of  t' 
on  the  right-hand  sides  of  these  equations  equal  to  1.  The  results  of 
these  integrations  are 

u'  =  |(l  -  r'2)  (44) 


and 
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(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


RESULTS  AND  DISCUSSION 
Laminar  Flow  of  Gases 

Velocity  and  temperature  distributions .  -  Velocity  and  temperature 
distributions  for  fully  developed  laminar  flow  of  gases  in  tubes  with 
variable  fluid  properties  are  shown  in  figures  1  and  2,  respectively. 

The  profiles  were  calculated  from  equations  (42)  and  (43),  where  the 
various  quantities  in  the  equations  were  calculated  as  indicated  after 
the  equations.  The  exponent  d  was  found  from  viscosity  data  to  have 
an  average  value  of  0.68  for  air  and  most  common  gases  for  temperatures 
between  0°  and  2000°  F . 

It  is  seen  from  figure  1  that  heat  addition  to  the  gas  sharpened 
the  peak  of  the  velocity  profile  in  the  central  portion  of  the  tube, 
whereas  heat  extraction  causes  a  flattening  of  the  profile  at  the  center 
of  the  tube.  These  results  agree  qualitatively  with  the  sketch  given  in 
reference  6  (p.186).  Inasmuch  as  the  viscosity  of  a  gas  increases  with  tem¬ 
perature,  the  viscosity  near  the  wall  for  heat  addition  is  greater  than 
that  near  the  center  of  the  tube  so  that  the  gas  is  slowed  down  near  the 
wall  and  accelerated  in  the  center  of  the  tuhe  for  heat  addition  to  the 
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gas.  The  opposite  effect  is  produced  for  heat  extraction  from  the  gas. 
The  variation  of  density  across  the  tube  apparently  has  a  smaller  effect 
on  the  velocity  distributions  than  does  the  variation  of  viscosity. 

The  changes  in  sharpness  of  the  temperature  profile  at  the  tube 
center  with  heat  addition  and  extraction,  as  indicated  in  figure  2,  are 
in  the  same  direction  as  those  for  velocity  profiles  because  the  thermal 
conductivity  varies  in  the  same  vay  •with  temperature  as  does  the  viscos¬ 
ity.  The  shapes  of  velocity  and  temperature  profiles  for  laminar  flow, 
however,  differ  somewhat  in  contrast  to  turbulent  flow  where  velocity 
and  temperature  profiles  are  similar  (reference  7).  Although  equa¬ 
tions  (l)  and  (16)  have  the  same  form,  the  shear  stress  varies  linearly 
with  radius  whereas  the  heat  transfer  per  unit  area  does  not,  and  hence 
the  shape  of  the  velocity  and  temperature  profiles  may  be  expected  to 
differ.  In  turbulent  flow  the  profiles  are  so  flat  that  It  makes  little 
difference  what  assumptions  are  used  for  the  variations  of  shear  stress 
and  heab  transfer  with  radius. 

Nusselt  numbers  and  friction  parameters .  -  Nusselt  numbers  and 
friction  parameters  are  plotted  as  functions  of  ‘bg/^b  fisu**68  3  and 
4,  respectively.  The  Nusselt  numbers  were  calculated  from  equa¬ 
tions  (32),  (33),  and  (36)  and  the  friction  parameters  were  calculated 
from  equations  (38)  and  (39) .  It  is  seen  from  the  figures  that  the 
parameters  vary  considerably  with  to/t^  when  the  fluid  properties  are 
evaluated  at  either  the  wall  temperature  or  the  bulk  temperature.  The 
effect  of  to/t-fc  on  the  parameters  can,  however,  be  nearly  eliminated 

by  evaluating  the  conductivity  in  the  Nusselt  number  at 
t_0  =  -0.27 (tQ-tfc)  +  %  and.  by  evaluating  the  viscosity  in  the  fric¬ 

tion  parameter  at  ^0.58  =  Q-58(tQ-t^)  +  t^.  The  fluid  properties  in 

the  Nusselt  number  must  be  evaluated  at  a  different  temperature  than  that 
at  which  the  fluid  properties  in  the  friction  parameter  are  evaluated, 
inasmuch  as  the  velocity  and  temperature  profiles  are  dissimilar  and  the 
bulk  velocity  and  the  bulk  temperature  are  weighted  differently  (see 
equations  (22)  and  (24)). 


Laminar  Flow  of  Liquid  Metals 

Velocity  and  temperature  distributions.  -  The  velocity  and  tempera¬ 
ture  distributions  for  laminar  flow  of  liquid  metals  shown  in  figures  5 
and  6  were  calculated  in  the  same  way  as  were  those  for  gases.  The 
exponent  d  was  found  from  the  limited  viscosity  data  available  for 
liquid  metals  to  have  an  average  value  of  about  -1.6  for  lead,  bismuth, 
and  sodium  potassium  alloy  at  temperatures  between  300°  and  1500°  1 . 

The  changes  in  velocity  -prof ile  shapes  shown  in  figure  5  are 
opposite  to  those  for  laminar  flow  of  gases,  because  the  variation 
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with,  temperature  of  the  viscosity  of  liquid  metals  is  opposite  to  that 
of  gases.  It  is  seen  in  figure  6  that  the  effect  of  tp/tb  on  the 
shape  of  the  temperature  profile  is  slight  for  liquid  metals,  because 
the  thermal  conductivity  is  considered  constant,  so  that  any  change  in 
the  temperature  distribution  must  be  caused  by  variation  of  viscosity. 

Nusselt  numbers  and  friction  parameters.  -  Nusselt  numbers  and 
friction  parameters  are  plotted  as  functions  of  tp/t^  in  figures  7  and 

8,  respectively. 

The  curves  in  figure  7  show  that  the  friction  parameter  can  be  made 
substantially  independent  of  tp/tb  by  evaluating  the  viscosity  at 
t0>54  =  0.54(t0-t-b)  +  tb.  This  temperature  is  close  to  the  temperature 
used  for  evaluating  the  viscosity  in  the  friction  parameters  for  gases. 

The  curve  in  figure  8  shows  the  variation  of  Nusselt  number  with 
tp/tb  for  liquid  metals.  The  effect  of  tp/tb  on  the  Nusselt  numbers 
cannot  be  eliminated  as  in  the  case  of  friction  parameters,  inasmuch  as 
the  thermal  conductivity  of  liquid  metals  is  assumed  not  to  vary  with 
temperature.  The  variation  of  Nusselt  number  with  tp/tb  is  caused  by 
the  variation  of  viscosity  with  temperature  and  can  be  nearly  eliminated 
by  multiplying  the  Nusselt  numbers  in  figure  8  by  the  empirical  correction 

for  liquids  usually  given  in  the  literature  ( lip/p.-^)^' (See,  for 

example,  reference  6,  p.  190.)  The  fact  that  the  Nusselt  number  for 
liquid  metals  varies  with  tp/tb  indicates  that  Nusselt  numbers  for  gases 

are  probably  affected  by  variation  of  viscosity  as  well  as  variation  of 
conductivity.  Evaluating  the  conductivity  at  a.  certain  temperature  in 
order  to  eliminate  the  effect  of  tQ/tb  seems  to  be  only  a  convenient  way 
of  presenting  the  results. 

Comparison  of  Laminar  and  Turbulent  Flow  of  Gases 

A  comparison  of  laminar  and  turbulent  flow  of  gases  is  given  in 
figures  9  and  10. 

In  figure  9  are  shown  calculated  laminar  Nusselt  numbers  from  the 
present  investigation  and  calculated  turbulent  Nusselt  numbers  for  a 
Prandtl  number  of  1  from  reference  7.  The  laminar  Nusselt  number  is 
independent  of  the  Reynolds  and  Prandtl  numbers,  whereas  the  turbulent 
Nusselt  number  is  dependent  on  them.  In  order  to  eliminate  the  effects 
of  tQ/tb,  the  fluid  properties  including  the  density  in  the  Nusselt  and 

Reynolds  numbers  are  evaluated  at  the  average  of  the  wall  and  bulk  temper¬ 
atures  tp.5  for  turbulent  heat  transfer  and  at  a  temperature  between 
the  bulk  temperature  and  the  temperature  at  the  center  of  the  tube  t_p.27 
for  laminar  heat  transfer.  The  fluid  properties  for  heat  transfer  in  the 
transition  region  should  evidently  be  evaluated  at  temperatures  between 
these  two. 
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In  figure  10,  the  laminar  friction  factors  were  calculated  from 
equation  (41)  and  the  turbulent  friction  factors  were  taken  from  refer¬ 
ence  7 .  The  laminar  friction  factor  f  is  independent  of  Prandtl 
number.  (The  laminar  friction  parameter  5  is  independent  of  both 
Reynolds  and  Prandtl  numbers . )  For  friction  factors ,  the  fluid  prop¬ 
erties  including  the  density  for  both  laminar  and  turbulent  flow  are 
evaluated  at  temperatures  close  to  the  average  of  the  wall  and  bulk 
temperatures . 


SUMMARY  OF  RESULTS 

The  following  results  were  obtained  from  an  analytical  investi¬ 
gation  of  laminar  flow  in  tubes : 

1.  Variation  of  fluid  properties  across  the  tubes  for  laminar  flow 
of  gases  caused  a  sharpening  of  the  peak  of  the  velocity  and  temperature 
profiles  at  the  center  of  the  tube  for  heat  addition  to  the  gas  and 

a  flattening  of  these  profiles  for  heat  extraction. 

2.  Variation  of  fluid  properties  across  the  tube  for  laminar  flow 
of  liquid  metals  caused  a  flattening  of  the  velocity  profile  at  the  cen¬ 
ter  of  the  tube  for  heat  addition  to  the  liquid  metal  and  a  sharpening 
of  the  profile  for  heat  extraction.  Variation  of  fluid  properties  had 
but  a  slight  effect  on  the  temperature  profile  for  liquid  metals. 

3.  The  effects  of  ratio  of  wall  to  bulk  temperature  on  This se It 
numbers  and  friction  parameters  for  laminar  flow  with  variable  fluid 
properties  could  be  eliminated  by  evaluating  the  fluid  properties  in 
these  quantities  at  various  temperatures  in  the  fluid. 

4.  The  velocity  and  temperature  profiles,  as  well  as  Nusselt  number 
and  friction  parameter,  were  independent  of  Reynolds  and  Prandtl  number 
for  fully  developed  laminar  flow  with  fluid  properties  that  varied 
across  the  tube. 


National  Advisory  Committee  for  Aeronautics, 
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Figure  1.  -  Velocity  distributions  for  laminar 
flow  of  gases  in  tubes. 


3.  -  Variation  of  Nusselt  number  with  ratio  of  wall  to  bulk  temperature  for  laminar  flow  of  gases 
with  conductivity  evaluated  at  temperatures  of  t_Q#27  =  -0.27  (tQ-t^)  +  t^,  t^,  and 


Figure  4.  -  Variation  of  friction  parameter  with  ratio  of  wall  to  bulk  temperature 
gases  with  viscosity  evaluated  at  temperatures  of  t^>  ^0.58  =  0.58  (tQ-t^) 
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Figure  6.  -  Temperature  distributions  for  laminar 
flow  of  liquid  metals  in  tubes. 
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Figure  7.  -  Variation  of  friction  parameter  with  ratio  of  wall  to  bulk  temperature  for  laminar 
medals  with  viscosity  evaluated  at  temperatures  of  to,  to  S4.  =  0.54  (to-tv)  +  tK.  ai 


metals  with  constant  thermal  conductivity. 
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Figure  9.  -  Variation  of  Nusselt  number  with  Reynolds  number  for  flow  of  gases  with 
fluid  properties  including  density  in  Nusselt  number  and  Reynolds  number  evaluated 
at  temperatures  of  t_oa27  “  -0*27  (tQ-t^)  +  for  laminar  heat  transfer  and  at 
temperatures  of  t0#5  =  0.5  (tp-t^)  +  tb  for  turbulent  heat  transfer  at  Prandtl 
number  of  1. 


